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Abstract 

We explore the dynamical stabilities of a quasi-one dimensional (ID) Bose-Einstein 
condensate (BEC) consisting of fixed N atoms with time-independent external potential. 
For the stationary states with zero flow density the general solution of the perturbed 
time evolution equation is constructed, and the stability criterions concerning the ini- 
tial conditions and system parameters are established. Taking the lattice potential case 
as an example, the stability and instability regions on the parameter space are found. 
The results suggest a method for selecting experimental parameters and adjusting initial 
conditions to suppress the instabilities. 
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1 Introduction 

Experimental observation of atomic gas Bose-Einstein condensates (BECs) has caused signif- 
icant stimulation to the study of macroscopic quantum phenomena with nonlinearity. In the 
* Author to whom all correspondence should be addressed, Email address: whhai2005@yahoo.com.cn 
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mean field regime where the BECs are governed by the Gross-Pitaevskii equations (GPE), the 
BEC of a stationary state can be observed carefully in experiments only for the stable solutions 
of GPE. For the purpose of applications, the studies on the stability and instability of the 
solutions of GPE are necessary and important pQ-pT]. Recently, the instabilities of BECs have 
attracted much interest and the corresponding experimental [HOE] and theoretical [i]-|13j 
works were reported for various BEC systems. Several different definitions such as the Landau 
instability dynamical instability [11], quantum instability [12], parametric instability [9] 

and modulational instability [18] were employed. The used research methods concerned the 
characteristic exponent technique pE], Gaussian variational approach [7], and the numerical 
simulations to the partial differential equations [5J [HI [TO]. The reported results showed that 
the instabilities are associated with the BEC collapse [TBI [19] , implosion and chaos [20] - [26J, 
dynamical superfluid-insulator transition [27] , and the formation and evolution of the matter- 
wave bright solitons [2H1 ESI ED]- In order to stabilize the BECs [17], some stability criteria 
[T5] and parameter regions [H EJ EH E2] were demonstrated. Most of the works focus in the 
stabilities under random perturbations. Experimentally [3| and theoretically [5] investigating 
the stabilities under the controllable perturbations has also become a challenging problem. 

In the sense of Lyapunov, the instability entails that the initially small deviations from the 
unperturbed state grow without upper limit. We shall restrict the dynamical instability to the 
particular case of nonzero characteristic exponents such that the minor deviations from the 
unperturbed state grow exponentially fast [U [TT] . All of the above-mentioned investigations 
on the dynamical stabilities and instabilities are based on such a type of instability. By the 
control of instability we mean to induce the transitions from unstable states to stable ones. 
Realization of the control needs selecting the system parameters to enter the stability regions, 
or initially using a controllable perturbation as a control signal to suppress the growth of 
perturbed solutions. Any experiment always contains a degree of noise, that leads to the 
random perturbations to the system. Therefore, in order to suppress the known unstable 
motions, we have to initially adjust the system by using the control signal being stronger than 
the noise. 

In the previous work, we have investigated the stabilities of BECs for the time-dependent 
chaotic states [221 ES] and dissipative cases [3TJ • In this paper, we shall consider the dynamical 
stability of the stationary states for a quasi- ID BEC consisting of fixed iV atoms with time- 
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independent external potential and atomic scattering length. It will be demonstrated that 
for the case of zero flow density the bounded perturbed solutions depend on the external 
potential, condensed atom number, and the initial disturbances. The dependence implies that 
the stationary state of BEC is certainly stable only for the given parameter region and the 
possible instability can be suppressed by some initial adjustments. We take the BECs held 
in an optical lattice as an exemplification to illustrate the results on the stability, instability 
and undetermined stability. The results contain the known analytical assertions for the optical 
potential case [5j[TT] and supply a method for selecting experimental parameters and adjusting 
initial conditions to establish the stable motions of BEC. 



2 Linearized equations and their solutions in the case of 
zero flow density 

We start with the dimensionless quasi- ID GPE [TH [33j EI] 

i^t = ~^xx+\y(x)+g 1 \i/)\ 2 ]i/) ) (1) 

where the suitable units with % — m — 1 have been considered, V(x) denotes the external 
potential, the quasi-lD interaction intensity g\ is related to the s-wave scattering length a s , 
atomic mass m and the transverse trap frequency uj r [351 136] for the normalized wave-function ijj 
with norm \ip\ 2 being the linear density of atomic number [TTII34"] . It is well known that different 
solutions of a nonlinear equation may possess different stabilities. Here we study stability only 
for the stationary state solution of the form 

■ipo = R(x) exp[i9(x) — i/it], (2) 

where is the chemical potential, R(x) and 9{x) represent the module and phase, which are 
both the real functions. In the considered units, the phase gradient 9 X is equal to the flow 
velocity field. Given the module, we define the useful Hermitian operators [UJ 

1 d 2 

L n = _ 2£fo2 + n 9iR 2 + V(x) - /j,, for n= 1,3. (3) 
Then inserting Eq. (2) into Eq. (1) gives the equations 

Li#(x) = 0. (4) 
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In the equation we have assumed the flow velocity field and current density being zero. 

We now investigate the stability of stationary state Eq. (2) by using the linear stability 
analysis, which is associated with boundedness of the perturbed solution [HI [T5] 

ip = [R(x) + efa(x, t) + ie0 2 (x, t)} exp[i9(x) — i/xt], (5) 

where the perturbed correction e4>i(x,t) is real function with constant |e| <C 1. Substituting 
Eqs. (5) and (4) into Eq. (1) yields the linearized equations [11] 
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Li4> 2 , 4>2t = ~L 3 fa. (6) 



For most of external potentials V(x) we cannot derive the exact solutions from Eq. (1) or 
Eq. (4) such that the operators L n cannot be determined exactly. In the case of optic lattice 
potential, some specially exact solutions have been found [HJ QUI [SB], however, solving Eq. 
(6) for the general solution is still difficult. Therefore, we have to focus our attentions to the 
dynamical stability which is associated with the perturbed solutions of space-time separation, 

<f>i{x,t) = Ti(t)ipi(x), for z = l,2. (7) 

Note that the real function fa limits and ipi to real or imaginary simultaneously, the difference 
between both is only a sign " — " of fa. We take real Ti, (pi, T 2 , and y? 2 without loss of generality, 
since the changes of the signs of fa do not affect the stability analysis. We shall discuss how to 
establish the sufficient conditions of stability as follows. 

Combining Eq. (6) with Eq. (7), we get the coupled ordinary differential equations 

W) = AxT 2 (t), f 2 {t) = -\ 2 T 1 {t)- (8) 
L 3 <Pi(x) = A 2 y3 2 (x), Lx<f 2 (x) = Xi<fx{x). (9) 

Here Aj is the real eigenvalue determined by the initial perturbations Tj(0), Tj(0). The corre- 
sponding decoupled equations are derived easily from the coupled ones as 

r,M= - a.a.t.w, A^lg, a, = -^ ; (io) 

L-iLz<f\ = AiA 2 y2i, L^L^ = AiA 2 ^ 2 . (11) 
Obviously, the general solutions of Eq. (10) can be written as the exponential functions 



Ti = Ae xt + B ie - Xt , A = v/-AiA 2 



At 



1^ 



A" 
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Ti(0) + -Ti(0) , Bi = - Ti(0) - -Ti(0) , (12) 



A" 



where A iy Bi are real or complex constants, which make Tj(t) the real functions. Based on the 
existence of bounded eigenstates tpi(x), the results are classified as the three cases: 

(j) Stability criterion: The eigenstates of Eq. (11) are bounded if and only if their 
eigenvalues are positive, A1A2 = —A 2 > 0, that makes A the imaginary constant and Tj the 
periodic functions. 

(jj) Instability criterion: One can find a negative eigenvalue AiA 2 = —A 2 < associated 
with a set of bounded eigenstates of Eq. (11) that makes Tj the real exponential function. 

(jjj) Undetermined stability: One cannot determine whether all eigenvalues of the bounded 
eigenstates of Eq. (11) are positive. In this case, we can use criterion (j) to control the possible 
instability of case (jj). 

3 Stability regions on the parameter space and control 
of instability 

It is interesting noting that if the initial perturbations can be determined, the dynamical 
instability of real A case can be controlled by adjusting the initial disturbances to obey = 
that will suppress the exponentially rapid growth of T in Eq. (12). From Eqs. (12) and (10) 
we establish the controlling criteria for the instability as T(0) = — AT(0). However, for the 
random initial perturbations such a control is difficult to do, since we cannot determine the 
initial values T(0) and T(0). Therefore, in the case of random perturbation we are interested 
in determining the same eigenvalue A1A2 of operators L1L3 and L3L1, since the stability can be 
established if and only if the eigenvalue is positive such that A 2 = — A X A 2 < 0. Let a > a g and 
P > flg be the eigenvalues of operators L\ and L 3 , which are determined by the eigenequations 
L\u[x) = au(x), Lsv(x) = /3v(x) with u and v being their eigenfunctions, where a g and f} g 
express the corresponding ground state eigenvalues respectively. From Eq. (3) we know the 
relation L 3 = L\ +2giR 2 that means a g < (3 g for g 1 > and a g > f3 g for g 1 < 0. It is clear that 
Eq. (4) is one of the eigenequations of L\ for the eigenvalue a = so that the ground state 
eigenvalue obeys a g < for any g±. Then f3 g can be positive or negative for g± > and j3 g < 
for g 1 < 0. 

From the above-mentioned results we establish the stability and instability conditions: 
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Case gi > 0: The sufficient condition of stability is a g = 0, since such a ground state 
eigenvalue implies a > and (3 > for all of the eigenstates such that the well known spectral 
theorem gives (HKTUIEZ] A1A2 > 0. The corresponding sufficient conditions of instability reads 
a g < and (3 g > 0. 

Case g\ < 0: The ground state eigenvalues satisfy the inequality (3 g < a g < 0. So the 
sufficient condition of instability is a g = 0. 

In all the other cases, we don't know whether A1A2 is certainly positive or negative, so the 
linear stabilities are analytically undetermined. It is worth noting that Eq. (4) infers R(x) to 
be one of the eigenstates of L\ with eigenvalue cur = 0. Therefore, if R(x) is a ground state, 
the above stability and instability conditions indicate that this state is stable for gi > 0, and 
unstable (or metastable) for g\ < 0. 

Note that all the above-mentioned results are valid for arbitrary time-independent 
potential. We will take the BEC held in an optical lattice as a concrete physical example to 
evidence these results. In the lattice potential case, the above-mentioned sufficient con- 
ditions agree with the stability and instability criterions established by the authors of Ref. 
[TTJ. We shall apply the sufficient stability and instability conditions to find the corresponding 
stability and instability regions on the parameter space, and apply these results to study the 
stabilization of the considered BEC system. 

For an arbitrary time-independent potential, the eigenequation L\u = au can be rewritten 
as the integral form [58] 




where q > is a real constant. This integral equation can be directly proved by taking the 
second derivative from its both sides. The integrals in Eq. (13) are indefinite, what means 
that the solutions are defined with accuracy of two additive constants. While the eigenequa- 
tion L\u = au just is a second order equation which also implies two arbitrary constants 
determined by the boundary conditions. It is the two additive constants to make the integral 
equation (13) completely equivalent to the eigenequation. The stability requires the eigenstate 
to be bounded and the possible bounded solution u must satisfy the boundedness condition 




u 




(13) 
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lim^-too J e Tqx fudx = 0. Under this condition and for the lattice potential case,we can 
apply the l'Hopital rule to get the superior limit [25j 

lim u < lim u x + lim w 2 = 2q~ 2 lim (fu). (14) 

Note that there is not the usual limit, because of the periodicity of lattice potential. It is clear 
that the solution of linear equation L\u = au can be taken as u(x) = Au'{x) with arbitrary 
constant A and any solution u'(x) such that one can always select u to obey \im. x ^± 00 u > 0. 
Thus Eq. (14) implies 2g -2 lmx r ^ = |- 00 / > 1, namely 

a>-{fM + ^^~[-V(x)-g 1 R 2 (x))}^a g . (15) 

For the eigenequation L 3 v = /3v after using 3g x instead of g x , the same calculations give 

P>-{fA + Jjjr[-V(x) - 3 9l R 2 (x)]} = /3 g . (16) 

Combining Eq. (15) with the stability sufficient condition [TTJ a g = for gi > 0, we get the 
parameter region of stability 

fi = H a = -Tim" [-V(x) - giR 2 (x)] for 9l > 0, (17) 

x — >±oo 

which contains the relation among fi,gi and the potential parameters. Applying Eqs. (15) and 
(16) to the instability sufficient conditions [11] a g < 0, (3 g > for g x > and a g = for g x < 0, 
we get the parameter regions of instability 

- lim [-V(x) - giR 2 (x)} < \i = \i in < - lim [-V(x) - 3g 1 R 2 (x)} forg 1 >0; 

n±oo x— >±oo 

fi in = -TiS~ [-V(x) - gi R 2 (x)} for 9l <0. (18) 

x — >±oo 

By the sufficient conditions we mean that the stationary state R(x)e~ lflt of Eq. (1) is certainly 
stable for the \i values in the parameter region fixed by Eq. (17), and the stationary states 
are certainly unstable for the \x values in any region of Eq. (18). The dynamical stabilities are 
undetermined outside. 

We now see the physical meaning of the stability relation in Eq. (17) for the stationary states 
of BEC with zero current density. Setting the sum of external potential and internal interaction 
as U(x) = V(x) + giR 2 (x) with periodic V(x) and bounded R(x), when U(x) > B is satisfied 
for all x values and a fixed constant B, Eq. (17) implies /j, s = B < U(x). Namely the sufficient 
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stability condition means that if the chemical potential is equal to the minimum of U(x), the 
considered states are certainly stable. For a known state the stability can be easily examined by 
using Eq. (17). We have tested the exact solutions given in Ref. [TT] for the potential V(x) = 
— Vosn 2 (x, k) and found that some of them have the instabilities and undetermined stabilities, 
where | Vol is the potential depth and sn(x, k) the Jacobian elliptic sine function with k being 
the modulus. Substituting one of the exact solutions, giR 2 (x) = —(1 + Vo//c 2 )[l — k 2 sn 2 (x, k)} 
with the potential depth — Vo > k 2 and chemical potential fi = —1 — Vo/k 2 + k 2 /2 [see Eq. 
(12) of Phys. Rev. E63, 036612(2001)], into Eq. (17) yields the stability parameter relation 
ix s = — 1 — Vo/k 2 . A difference of k 2 /2 exists between the fi s value required by the stability 
condition and the chemical potential \x in the exact solution, namely the stability criterion (17) 
is not met here. This assertion differs from the result of Ref. [llj . where this solution fits their 
stability criterion and the stability is independent of the parameters k and Vo- However, when 
the potential depth |Vo| is much greater than the modulus k (e.g. Vo = —1 and k = 0.2), we 
have the chemical potential near the stability relation (17) (// = 24.02 = fi s + 0.02 // s ). This 
infers the higher stability being associated with a smaller value of the modulus k and a relatively 
greater |Vo| value. Thus our stability parameter criterion suggests that for a known solution 
with instability or undetermined stability one can raise the practical stability by adjusting the 
system parameter (e.g. the above k and |Vo|) to approach the values of the stability region in 



Generally, constructing a stable exact solution of GPE is not easy, because of the non- 
integrability of Eq. (4) with periodic potential. However, in the large- N limit, we can fulfil the 
criterion (17) for the case of a repulsive nonlinearity, since the Thomas- Fermi (TF) approxi- 
mation [33J U(x) = fiTF just fits the stability relation. Therefore, it is practical relevant to 
prepare such a stable TF state R(x, Htf) by increasing the condensed atom number N. Given 
the number N and the periodic boundary condition experimentally, from the normalization 
condition N = n f£ R 2 (x, fiTF)dx = n Jq[^tf — V(x)]dx/gi we derive the chemical potential of 
the stable TF state 



which is related to the atom number N and the potential strength Vo and period K(k), where 
n ~ 100 is the lattice number. In fact, noticing the dependence of R = R(x, //) on \i in Eq. (4), 



Eq. (17). 




(19) 
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the normalization condition of any known state can also lead to \i — /i(A r ) and R = R(x,N). 
Applying them to eliminate /i in Eqs. (17) and (18) will give the corresponding relationships 
among the experimental parameters N, gi, Vq and K(k). So we can control the instability of 
the known state by selecting the experimental parameters to fit or to approach the stability 
region of Eq. (17). In many practical cases, we cannot obtain the exact solution of Eq. (4) 
for some periodic potentials, that necessitates the numerical investigation. In order to fit (or 
near) the stability region in Eq. (17) and to avoid the instability regions of Eq. (18), we could 
use Eq. (19) to estimate and adjusted the chemical potential in region /i ps \itf such that the 
stability of the numerical solutions of Eq. (4) can also be established or improved. 

On the other hand, in the case of arbitrary time-independent potential, for some known 
unstable solutions R = R(x,/ii n ) from Eqs. (10) and (12) we can experimentally set and adjust 
the initially controllable perturbation as a control signal [3J to suppress the exponentially fast 
growth of Ti(t). Although the phase 9 and amplitude R are time-independent in the considered 
case, the initial perturbations can result in the nontrivial and time-dependent corrections to 
the phase and atomic-number density. From Eq. (5) we find their first corrections as 

A9(x,t) « aTctan[ET 2 (t)ip 2 (x)/R(x)] w eT 2 (t)<p 2 (x)/R(x), 

A\iP\ 2 (x,t) « 2eT 1 (t)<p 1 (x)R(x), (20) 

which are initially proportional to Xi(0) and T 2 (0) respectively. Making use of Eq. (20), 
the adjustments to the initially controllable perturbations can be performed by trimming the 
number density \ip\ 2 , velocity field (A6) x and their time derivatives which are proportional to 
the corresponding trimming velocities. Given Eqs. (10) and (12), we know the stability initial 
criterion 

A 2 = -AiA 2 = 7\(0)T 2 (0)/pi(0)T 2 (0)] < 0. (21) 

Once Eq. (21) is satisfied in the adjustments to the initial perturbations, Eq. (12) becomes the 
periodic solution which implies the stability. Although we cannot determine the initial values 
Tj(0) and Tj(0), experimentally, the number density can be adjusted by varying the condensed 
atom number, and the adjustments to superfluid velocity may be related to a displacement 
Ax of a magnetic potential [3J. According to Eqs. (20) and (21), if we initially increases (or 
decreases) both the relative derivative T 2 (0)/T 2 (0) = dAe ^ x '^ \ t=0 / A8 x (x, 0) of flow velocity and 
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the relative derivative T 1 (0)/T 1 (0) = aA M fo*) | t=o /A|'0| 2 (a;, 0) of atomic number density, the 
stability initial criterion (21) is destroyed and the system will become unstable. But when one of 
them is increased and another is decreased simultaneously, the stability criterion (21) is satisfied 
and the possible instability is suppressed. These assertions may be tasted experimentally. 

4 Conclusions and discussions 

In conclusion, we have investigated the dynamical stability, instability and undetermined sta- 
bility of a quasi-lD BEC in the stationary states for time-independent external potential and 
atomic scattering length, and fixed atomic number. After space-time separation of variables, 
we derive the general solutions of the linearized time-evolution equations for the trivial phase 
case and give a stability criterion related to the initial conditions. As an important example, we 
evidence the stability criterion analytically for the BEC held in an optical lattice potential. By 
using the known sufficient conditions of stability and instability [TT], several parameter regions 
of stability and instability are shown. Our results contain some new stability predictions which 
can be tested with current experimental setups. Finally, we stress that applying our stability 
initial criterion and parameter region one can stabilize the considered BEC system by adjusting 
the system parameters experimentally to enter or near the stability region of Eq. (17) on the 
parameter space. For the parameters out of the stability region we can also establish or improve 
the stability by adjusting the initial flow velocity and atomic number density to fit or approach 
the stability initial criterion. 
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